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Using Liischer's method we determine the elastic scattering phases in the broken phase of the 4-dimensional 
0(4) non-linear cr-model from the two-particle energy spectrum in a Monte-Carlo study on finite lattices. In 
the isospin-O-channel we observe the cr-resonance and extract its mass and its width. In all scattering channels 
investigated the results are consistent with perturbative calculations. 
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1. INTRODUCTION 

Liischer established a relation between the energy 
spectrum of two-particle states in a finite box 
with periodic boundary conditions and elastic 
scattering phase shifts defined in infinite volume 
[1]. Since two-particle energy levels are calculable 
by Monte Carlo techniques, this relation opens 
the possibility to extract phase shifts from nu- 
merical simulations on finite lattices. We give 
a resume of our results on the elastic scattering 
phases in the broken phase of the 4-dimensional 
0(4) non- linear cr-model including the corrections 
to our preliminary results of the Lattice '92 con- 
tribution [2]. For a more detailed discussion see 
ref. [3]. 

2. FROM THE TWO-PARTICLE SPEC- 
TRUM TO THE SCATTERING PHA- 
SES 

We give the "pions" of our model a nonzero mass 
by means of an external source J in the ac- 
tion: 



5=-2«:^^$^$^+^ + J^$^ (1) 

The scalar field is represented as a four- 
component vector of unit length: ^^^^ = 1. 
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Two-pion states are classified according to cu- 
bic symmetry and have "isospin" 7 = 0, 1, 2. For 
their investigation wc define the operator: 



fi ,t 



(2) 
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where /^(ti), i — 1,2,... is some basis of wave 
functions with correct cubic symmetry and ^ 
is the spatial Fourier transform of the field $g ^ 
on a lattice of spatial extent L: 



3>1 
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(3) 



The simplest choice of cubically invariant wave 
functions is a sum over plane waves: 



fi(n) = 6i^ft2 



0,1,2,. 



(4) 



We also used another set of (Liischer-) wave func- 
tions [1], which gave the same final results (see [3] 
for further details). 

In the isospin-0 channel, where the cr-resonance 
is expected, we use operators Oi{t) given by 



1 ^ 



(5) 



Additionally we have to take into account the cr 
field at zero momentum 



(6) 
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since it has the correct quantum numbers and is 
expected to create a state with energy below the 
inelastic threshold. 

The two-particle energies W^, are extracted 
from the matrices of connected correlation func- 
tions [4]. For the different isospin channels these 
matrices are given by 

C^t) = (O,(i)O,(0))„ 

Cj.{t) = {lmOf(t)lmOf{0)),, (7) 

Cf^it) = (ReOf(i)Re0f(O)),-CO.(f). 

As described in [1], for each lattice extent L 
and each two-particle energy level < 4771^ we 
get one value of the scattering phase shifts 6q, 
S\, and ^g, respectively, in the elastic region < 
j 771^ < \f?>. It is computed from the key 
relation 

5/(M = -'^(^) modTT, (8) 

where ^ is a continuous function defined in [1], see 
also [2,3] (Z=angular momentum). The momen- 
tum fc„ corresponding to has to be calculated 
with the help of the (lattice) energy momentum 
relation 

(2sinh^^^|^)'=777^ + C (9) 

This relation does not determine uniquely, 

since fe„ depends not only on fe„ but also on the 
direction of fc^. We include this small ambiguity 
in the error estimate for fc^. 

3. NUMERICAL RESULTS 

The data points in figs. 1- 3 show the momenta ex- 
tracted from the energy levels and the correspon- 
ding phase shifts 6;S,^ii,^g at (k = 0.315, J = 0.01). 

How do our results compare with one-loop per- 
turbation theory if we insert the renormalized 
coupling constant and masses as determined by 
our simulations ? The dashed curves are the per- 
turbative predictions based on an estimate vn^ 
of the resonance mass 777^, which was obtained 
from a fit to the non-resonance cr-propagator in 
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Figure 1. Comparison of the perturbative fit 
3. A. (solid curve) and the perturbative prediction 
(dashed curve) in the isospin-0 channel. 



momentum space. We consider tti^ to be an esti- 
mate only, since the a particle is unstable. Fig. 1 
shows that indeed 771^ lies below tti^. 

In figs. 2, 3 we observe no significant deviation 
from perturbation theory: For isospin 1 and 2 
the perturbative predictions depend on 777^ only 
weakly. 

In order to determine the resonance mass 777^ 
and the decay width from the measured scatte- 
ring phases we have employed different methods: 

A. Perturbative Fit 

The one-loop perturbative formula for Sq [5] (see 
[3] for details) as function of fc; tti^, 777^, can 
be used as a fit ansatz for all points in the elastic 
region. This leads to 777^ = 0.691(3) and = 
0.112(7) in good agreement with the results from 
the Breit-Wigner fit below. 

B. Breit-Wigner Fit 

Fitting to the (relativistic) Breit-Wigner formula 

has the major advantage of being free of addi- 
tional asumptions, but can only be applied near 
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Figure 2. Is there a p in the isospin-1 channel? - 
The scattering phase shifts for isospin 1 indicate 
that there is none as expected from perturbation 
theory. 



the resonance at 5q = 7r/2, so only few data 
points are used. Nevertheless we get for the reso- 
nance mass rUa and width Va: 771^ = 0.706(2) and 
= 0.130(9). These numbers have to be compa- 
red with the estimate ifia =0.720(1) and the per- 
turbative prediction rCT= Tjjha, . . .) =0.121(1). 

4. DISCUSSION AND CONCLUSIONS 

The results for the decay width agree reasonably 
well with the perturbative predictions, while the 
resonance masses vfia lie systematically (about 
5 %) below the estimates fh„ from the fit to the 
propagator in momentum space. However, this 
discrepancy should not be too surprising, because 
the width of the cr-resonance for our choice of pa- 
rameters is rather large: « 0.15 

Once again, renormalized perturbation theory 
has turned out to be very reliable in the four- 
dimensional ^"'-theory. Furthermore, we have de- 
monstrated the applicability of Liischer's method 
for studying particle scattering processes in mas- 
sive quantum field theories on finite lattices - at 
least for this model. 
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Figure 3. Scattering phase shifts for isospin 2. 
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